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We investigated the electronic and magnetic structure of Sr3_j;Caa;Ru207 (0 < x < 3) on 
the basis of the double-layered three-dimensional multiband Hubbard model with spin-orbit 
interaction. In our model, lattice distortion is implemented as the modulation of transfer 
integrals or a crystal field. The most stable states are estimated within the unrestricted 
Hartree-Fock approximation, in which the colinear spin configurations with five different 
spin-quantization axes are adopted as candidates. The obtained spin structures for some 
particular lattice distortions are consistent with the neutron diffraction results for Ca3Ru207. 
Also, some magnetic phase transitions can occur due to changes in lattice distortion. These 
results facilitate the comprehensive understanding of the phase diagram of Sr3_:rCa2.Ru207. 

KEYWORDS: three-dimensional multiband Hubbard model, double-layered ruthenate, lattice 
distortion, unrestricted Hartree-Fock, spin-orbit interaction 

1. Introduction 

The series of double- layered ruthenates, Sr3_a;Caa;Ru207 {0 < x < 3), possesses a va- 
riety of phases under a magnetic field (H) or pressure (P). One end- member of the series, 
Sr3Ru207, shows the magnetic field-tuned quantum criticality, which is accompanied by the 
metamagnetic transition around H ~ 7.85T.^) The metamagnetic transition, which was ini- 
tiahy observed by Cao et al.,^^ was revealed to be a double transition, and its second transition 
in the higher- field is sensitive to the field angle. Without an applied magnetic field, Sr3Ru207 
behaves as a Fermi liquid at low temperature.^^ Angle-resolved photo-emission spectroscopy 
(ARPES) has been used to observe its well-defined Fermi surfaces, and neutron diffraction 
methods revealed a lack of that long-range magnetic order. ^-^ Meanwhile, inelastic neutron 
scattering has been used to observe two-dimensional ferromagnetic fluctuations as incommen- 
surate peaks attributed to Fermi surface nesting.''^ These fluctuations induce ferromagnetism 
when uniaxial pressures along the c-axis are applied.^' Since neutron diffraction analysis 
displays the temperature and pressure effects on the crystal structure, the phase transition 
induced by uniaxial pressures suggests that the ferromagnetic fluctuations are susceptible to 
structural changes. This situation is similar to the ferromagnetic ground state at the surface 
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of Sr2Ru04, as observed by scanning tunneling microscopy (STM).^i) In Sr2Ru04, the ferro- 
magnetic ground state arises due to a perturbative lattice distortion at the surface, that is, 
an in-plane rotation of the RuOg octahedron produced by the surface strain. 

Another end-member of the series, Ca3Ru207, shows two transitions at Tm = 48K and 
T/v = 56K.^^^ While Tjv has been confirmed as an antifcrromagnctic ordering temperature,^^) 
Tm was believed to be a Mott-likc metal-insulator transition temperature though quantum 
oscillations in the c-axis resistivity pc, for H \\ c are observed below Tm-^^'^^^ The electri- 
cal resistivity and optical conductivity spectra of single crystals grown by a floating-zone 
method proved that the ground state of Ca3Ru207 is quasi-two-dimensional metallic. ^^'^'''^ 
Furthermore, the magnetostriction data for the single crystals demonstrated that the first- 
order transition at Tm can be attributed to a discontinuous change in the lattice constants. 
The quantum oscillations observed in Ca3Ru207 have also shown that its ground state is 
metallic with low-carrier density. The magnetic structure of Ca3Ru207 in the ground state 
was clarified using neutron diffraction analysis: the magnetic moments align ferromagnetically 
within the double layer and antiferromagnetically between the double layers. While these 
magnetic moments lie along the 6- axis for T < Tm, the first-order transition at Tm changes 
their directions to align with the a-axis for T > Tm-"^^'"^"^^ Moreover, the two transitions at T/v 
and Tm are respectively weakened by the pressures along the c-axis and those within the aft- 
plane. These results indicate that the magnetic properties of Ca3Ru207 are also susceptible 
to structural changes. 

The concentration range < x < 3 of the Sr3_xCa'xRu207 series has been investigated in 
addition to the end members.^*^"^"-* For the intermediate x of this range, the system exhibits 
a variety of spin structures: a cluster spin-glass phase for 0.24 ^ x < 1.2,^^ •^'^^ and a canted 
antiferromagnetic phase for 1.2 < x < 2.0.^'"''' By comparing the Weiss temperatures for H \\ ah 
with those for -ff || c, Iwata et al.^^^ elucidate that the magnetic easy axis changes continuously 
from the a6-plane to the c-axis with decreasing x. Peng et al.'^^-* confirm that the magnetic 
easy axis is the 6-axis for x = 2.4 and x = 3.0. This result for x = 3.0 is consistent with the 
result of a neutron diffraction analysis for Ca3Ru207.^*^) Moreover, it has been reported that 
lattice constants vary with X in Sr3_^Ca^Ru207.2^'^°) 

Sr3_a;Caa;Ru207 (0 < X < 3) has also attracted a great deal of theoretical interest. The 
band structures of its end-members, Sr3Ru207 and Ca3Ru207, have been investigated with 
a local density approximation^^) or with the local spin density approximation (LSDA). 32,33) 
In particular, the magnetic field-tuned metamagnetic transition of Sr3Ru207 has been inten- 
sively studied as one of the electronic nematic phase transitions on the basis of microscopic 
theories. 3^"^°) The field-induced orbital-ordered phase of Ca3Ru207 was also investigated 
on the basis of the spin/orbital model.*^) However, few theoretical studies have investigated 
Sr3_a;Caa;Ru207 (0 < a; < 3), while a number of theoretical analyses have been performed 
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for the series of single-layered ruthenates, ^^-49) Ca2-a;Sr2:Ru04 (0 < a; < 2), which exhibit 
the Mott transition at a; ~ 0.2.^'') The two-dimensional multiband Hubbard model has been 
utilized by these theoretical analyses of the series of single-layered ruthenates. Meanwhile, in 
order to understand the series of double-layered ruthenates, we need to consider its three- 
dimensionality. The intrinsic importance of the three-dimensionality can be easily found in 
the experimental results introduced above (e.g., the magnetic structure of Ca3Ru207). In this 
paper, we investigate the electronic and magnetic structures of the double-layered ruthen- 
ate Sr3_a;Caa;Ru207 (0 < x < 3) on the basis of the three-dimensional (3D) multiband 
Hubbard model. Fully considering possible unequivalent sites and the spin-orbit interaction 
(SOI), we determine the ground state of our model for each lattice distortion within the 
unrestricted Hartree-Fock (UHF) approximation. Then, we find that the change in lattice 
distortion severely affects the electronic and magnetic structures. Our results suggest that the 
many physical phenomena of Sr3_a;Caa;Ru207 (0 < a; < 3) have a critical relationship with 
the change in the lattice distortion. 

2. Formulation 



Our 3D hubbard model with lattice distortion (Fig. 1) consists of A and B sub lattices. 
We consider every three t2g orbitals of the Ru 4d electrons located in these sublattices on 
the i-th layer (i = 1,2,3,4). Thus, our 3D Hubbard model Hamiltonian, H, is composed as 
follows: 



H 



1=1 j=l k 

4 

=1 k aa' 

4 



E E E l^lkJaa'^ika' + Bj^Jaa'^ika 
1=1 k aa' 

4 

- E E E [AWAfea + BiAk, 



(1) 



=1 k 



Here we use the abbreviations A, 



t 

iko 



.yz] .zx] .xy]'' 
^ika ^ika ^ika . 



A; = u Ay"" A'^y \ = 

J J , ^ika — y^ika ^ika ^ika) ' °ika — 

Bt^Bl^lB^^;)^ and B,k. - \^LBtL^L) > where ^,.(^1) and 5f,,(SfJj are the 
annihilation (creation) operators for the electron in the A and B sublattice on the i-th layer 
{i = 1,2,3,4), as specified by orbital ip = {yz, zx,xy}, momentum k, and spin a = {tji}, 
respectively, /x is the chemical potential. The nonvanishing h^^, hf^, hf^ki ^^fe ™ ®*1- (-*-) 
are 

/ Afc ^ 

Afe (z = l,2,3,4), (2) 
V efe / 



uAA _ iBB _ 
'Hik ~ 'Hik ~ 



3/16 



J. Phys. Soc. Jpn. Full Paper 




Fig. 1. (Color online) Unit cell of our 3D Hubbard model: (left) projection onto the a'c-plane, (center) 
projection onto the a'6'(a6)-plane, and (right) projection onto the 6'c-plane. The diamonds with 
solid circles(squares) represent the RuOe octahedrons on the A(B) sublattice. The Ru sites in 
these octahedrons are indicated by solid circles. Li indicates the i-th layer (j — 1, 2, 3, 4). and 6 
represent the rotation and tilting angles of the RuOg octahedron, respectively. 
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In eqs. (7)-(9) t± and t'j_ represent inter-layer transfers, and in eqs. (10)-(14), ti, t2, ts, t^, and 
Ao represent intra-layer transfers. A in eq. (13) represents the energy level difference between 
the dxy and dyz{dzx) orbitals due to the crystal field. The terms l^^, and l^^, in eq. (1), arising 
from the spin-orbit interaction, are determined from formulas that depend on the choice of 
the spin-quantization axis. Here we only consider coUinear spin states, with the five different 
spin-quantization axes as candidates for the most stable states. These five are the c-, a'-, h'-, 
a-, and 6-axes. When we consider the state with the spin-quantization axis parallel to the 
c-axis, we represent 1^^, and 1^^, by f^j:^^ and 1^^^'' , respectively. These are defined as follows: 
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Furthermore, when we consider the state with the spin-quantization axis parallel to the a-axis, 
we have 
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The interacting part H' in eq. (1) is represented by 
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where N is the number of fe-space points in the first Brillouin zone (FBZ). Here we only con- 
sider the on-site interactions, i.e. Coulomb repulsion in the same orbital U, Coulomb repulsion 
in different orbitals V, and the exchange interaction J. We adopt the UHF approximation 
with respect to every two sublattices, four layers, three orbitals and two spin states. Thus, we 
define 
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In order to determine the most stable of the five candidates, we conduct a self-consistent 
calculation for each candidate and find the one with the lowest electronic energy as estimated 
by eqs. (1) and (28). Hence, we can translate the most stable spin-quantization axis as the 
magnetic easy axis. Moreover, for the electronic state that wc determine has the most stable 
spin-quantization axis, we can calculate the five types of the magnetic order parameters. 
Four of these types are antiferromagnetic order parameters, referred to as Ai-AFM, A2-AFM, 
Ci-AFM, and C2-AFM and expressed as 
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h t2 h U t± t'j_ Xo A C U V J 

0.40 0.08 0.40 0.04 0.24 0.04 0.04 0.00 0.16 0.80 0.40 0.20 Fig. 2 

0.40 0.08 0.40 0.04 0.24 0.04 0.04 0.16 0.16 0.80 0.40 0.20 Fig. 3 

0.40 0.08 0.40 0.04 0.24 0.04 0.04 0.00 0.16 1.00 0.50 0.25 Fig. 4 

0.40 0.08 0.40 0.04 0.24 0.04 0.04 0.16 0.16 1.00 0.50 0.25 Fig. 5 



Table I. The parameter sets for the calculations. The parameter uint is eV. 
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respectively. The fifth type is the ferromagnetic order parameter, expressed as 

4 



m(FM) 



m. 



1=1 



Here, we introduce the magnetic momentum on each site: 
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Then we determne the magnetic order with the largest magnitude of these five order param- 
eters for each parameter set. 

3. Results and Discussion 

In our numerical calculations, we divide FBZ into a 20 x 20 x 20 equally spaced meshes. 
The parameter sets in eqs. (7)-(20) and (25) for these calculations are selected as shown in 
Table I. The choice of these parameters is based on preceding theoretical works. '^^^ The ratios 
V/U and J/U in each set are fixed at 0.5 and 0.25, respectively. For each set in Table I, both 
the rotation angle (p and the tilting angle 9 are varied from 0° to 20° by 5°. All resulting 
self-consistent fields and nf^^ have four digits of accuracy, and they satisfy 

EEE«<^ + <<^) = 32, (35) 

i=l ip a 

which means that there are four electrons per Ru site. 

We summarize our numerical results, as indicated by the last column of Table I, in Figs. 2- 
5. The figures show the magnetic easy axis, magnetic phase, and density of states (DOS's) 
determined for each pair of {(f), 6). Here, the magnetic phase is defined as the paramagnetism 
(PM) when none of the order parameters reach a finite value, and the magnetic phase is 
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Fig. 2. (Color online) The magnetic easy axis, magnetic phase, and density of states (DOS) deter- 
mined when A — 0.00 eV and U = 0.8 eV. The unit of each {(p, 9) is provided in "(degrees). Positive 
DOS is for spin up and negative DOS is for spin down. 



defined as X when several order parameters simultaneously reach the maximum value of the 
five. We can easily recognize that the various magnetic phases appear and that their phase 
transition is caused by the fractional change of the lattice distortion. The magnetic easy axis 
can vary even in the same magnetic phase. The variations in the magnetic phase and easy axis 
are caused by the transfers and SOI, which both depend on the lattice distortion. While the 
SOI dependence on the lattice distortion plays the primary role in the determination of the 
magnetic easy axis, the transfers dependence is more responsible for determining the magnetic 
phase. 

The energy level difference between dxy and dyz{dzx), i-e. A, also affects the electronic 
states since it relatively lowers the bands from the dxy orbital as well as enhances the full 
bandwidth W. When we compare Fig. 2 with Fig. 3, or Fig. 4 with Fig. 5, we find that a positive 
A makes the PM phase stable for a wider range of lattice distortion parameters. This can be 
derived from the decrease in U/W in conjunction with the change in A. In Ca3Ru207, the 
lattice constants abruptly change at Tm = 48K, where the first-order transition occurs"*^^' 
due to Jahn- Teller distortions of the RuOe octahedra.^^^ Thus, we can naturally assign our 
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Fig. 3. (Color online) The magnetic easy axis, magnetic phase, and DOS determined when A = 
0.16 eV and U = 0.8 eV. 



results for A > (Figs. 3 and 5) to the quasi-two-dimensional metallic state of Ca3Ru207 for 

In our model, all bands can be categorized into two types: bands derived from the dxy 
orbital or bands derived from the dyzid^x) orbital. When we respectively define their band- 
widths as Wxy and Wy^^zx, their dependence on (p and 9 is like Wxy oc (cos 20 cos and 
Wyz,zx Of (cos 0COS 2^)^, due to eqs. (10)-(12). When either the rotation {(p) or the tilting 
(6) increases, both Wxy and Wyz,zx decrease and U/W increases. Thus, the lattice distortion 
prefers the AFM phase rather than the PM phase due to a large U/W, as shown in Figs. 2 and 
3. Moreover, each bandwidth dependence creates differences in the lattice distortion effects 
on the electronic state between (p and 6. The results for ((/>, 6) = (20, 0) and (</>, 9) = (0, 20) in 
Figs. 2 and 3 provide clear evidence of these differences. 

Figs. 6 and 7 show the obtained antiferromagnetic structures. In these structures, Ai- 
AFM along the 6-axis is consistent with the structure of Ca3Ru207 observed by a neutron 
diffraction analysis. ^'^^ The Ai-AFM phase appears as shown in Figs. 2-3. This phase has been 
proved to be the most stable state by Singh and Auluck in the LSDA for Ca3Ru207, where 
it was noted as AFl.^^^ 
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Fig. 4. (Color online) The magnetic easy axis, magnetic phase, and DOS determined when A 
0.00 eV and U = 1.0 eV. 



The FM phase sometimes appears in the neighborhood of the Ai-AFM phase as shown in 
Figs. 2 and 4. In the Ai-AFM phase, the magnetic moments align ferromagnetically within 
the double layer, and the ferromagnetic correlation within this double layer is stronger than 
the antiferromagnetic correlation; this has been confirmed by an inelastic neutron scattering 
study for the spin- wave excitation in Ca3Ru207.^^^ In the intermediate regime where the 
AFM correlation between the double layer is not fully developed, electrons show the FM 
order as a whole although their magnetic moments are small. Thus, FM can be derived from 
the perturbative change of lattice distortion in our model, and this supports the emergence 
of FM induced by uniaxial pressures along the c-axis in Sr3Ru207.^'^) 

Let us note that the X magnetic phase appears in Figs. 4 and 5. Several spin configurations 
have almost the equivalent energy in this phase, so that the energy profile of this phase has 
a multi-valley structure in its ground state. Thus, we can identify the X magnetic phase as 
a cluster spin-glass phase in Sr3_2,Ca2,Ru207 for 0.24 ^ x < 1.2.^^"^^^ The variation of x in 
Sr3_^Caa;Ru207 is accompanied by the change in lattice distortion, ^'^''^^^ which must be one 
reason why the electronic state changes with x unless the carrier density remains the same. 
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Fig. 5. (Color online) The magnetic easy axis, magnetic phase, and DOS determined when A 
0.16 eV and U = 1.0 eV. 



4. Conclusion 

In this paper, we examined the lattice distortion effects on Sr3_i.CaKRu207 using the 
double-layered 3D multiband Hubbard model with SOI by the unrestricted Hartree-Fock 
calculation. For some types of lattice distortion, we obtained the Ai-AFM phase along the 
6-axis, consistent with the neutron scattering result for Ca3Ru207. Our results also indicate 
a possible ferromagnetic transition which is susceptible to the change in lattice distortion and 
the existence of a cluster spin-glass phase for the intermediate x. The electronic states with 
these above results are all metallic. This suggests that a number of physical phenomena in 
zero-field can be explained without the existence of a metal-insulator transition. In a recent 
experiment it was reported that the field-induced metamagnetic transition drives small lattice 
distortion in Sr3Ru207.^^^ To elucidate such a relation between lattice distortion and quantum 
critical phenomena, the study of the lattice distortion effects on Sr3_a;Caa;Ru207 in a finite 
magnetic field is needed in the future. 
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Fig. 6. (Color online) Antiferromagnetic structures obtained in our calculations. The black(red) solid 
circles represent the Ru site on the A(B) sublattice. 
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